June, 1994 



CMU-HEP94-20 
DOE/ER/40682-74 
PITT-94-06 



hep-ph/ 9406422 



Thermal Activation Rates in the Chirally Asymmetric 

Gross-Neveu Model 

Daniel BoyanovskyQ^ , David E. BrahmQ( 6) , 
Richard Holman^^, and Da-Shin LeefJ") 

^ Dept. of Physics and Astronomy, University of Pittsburgh, Pittsburgh PA 15260 
^Carnegie Mellon Physics Dept., Pittsburgh PA 15213 



Abstract 

We address the problem of how to incorporate quantum effects into the calculation 
of finite-temperature decay rates for a metastable state of a quantum field theory. To 
do this, we consider the Gross-Neveu model with an explicit chiral symmetry breaking 
term, which allows for a metastable state. This theory can be shown to have a "critical 
bubble" which is a solution to the exact equations of motions (i.e. to all orders in 
perturbation theory, including all higher derivative, quantum and thermal corrections). 
This configuration mediates the thermal activation of the metastable vacuum to the 
true ground state, with a decay rate V oc exp(— F c /T), where F c is the free energy of the 
critical bubble. We then compare this exact calculation to various approximations that 
have been used in previous work. We find that these approximations all overestimate 
the activation rate. Furthermore, we study the effect of finite baryon number upon 
the bubble profile and the activation barriers. We find that beyond a critical baryon 
number the activation barriers disappear altogether. 
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1 Introduction and Motivation 



Thermally activated reactions appear in many contexts, including chemical reactions [l|], cos- 
mological phase transitions and sphaleron-mediated baryon number violating transitions 
in the early universe^]. The physically important quantity in these transitions is the acti- 
vation rate; this sets the time scale for the process, and determines whether these reactions 
can be in local thermal equilibrium. 

There is a standard way of computing the activation rate in a field theory at finite 
temperature. Consider a field theory with a metastable false vacuum state as well as a 
stable ground state. The first step in the calculation is to construct a static, extremal 
configuration, corresponding to a critical bubble, in which the field is near the true vacuum 
on the inside, and in the false vacuum outside the bubble. Once this configuration is obtained, 
the activation rate T is proportional to exp(— F c /T), where T is the temperature, and F c is 
the free energy of the critical bubble. 

In most applications, the bubble configuration used in practice is a static extremum of the 
classical action^, However, this leaves open the question of how to incorporate quantum 
effects into the rate calculation ||. In particular, it is not clear whether the full theory has 
a critical bubble. Even if it does, the relation between the rate obtained by the using the 
classical bubble and the rate we would calculate using the exact bubble is not altogether 
clear, to say the least. 

In a general field theory, it is usually impossible to construct the exact critical bubble of 
the theory, so that the questions asked above are moot. Nonetheless, it would be of great 
interest to have some idea of how the exact and the approximate rates compare. It is for 
this reason that we consider the Gross- Neveu|7]] (GN) model in 1 + 1 dimensions. In the 
presence of a chiral symmetry breaking term, the theory has a metastable ground state. Of 
greater use to us, however, is the fact that we can obtain the exact critical bubble of the 
full quantum theory! Being able to do this allows us to proceed with the program described 
above. In particular, there are no "missing loops", or higher derivative corrections. While 
this is a theory in 1 + 1 dimensions, the fact that it has much of the same richness as 3 + 1 
theories (e.g. asymptotic freedom, dimensional transmutation, etc.) leads us to think that 
some of the results of this paper could be applicable to realistic field theories with metastable 
vacua in 3 + 1 dimensions. 

In the next section, we describe the Gross-Neveu model in detail. In particular, we 
exhibit the phase structure of the theory in the presence of an explicit chiral symmetry 
breaking term, and construct the exact bubble for this theory. We then proceed to compare 
the exponential part of the activation rate obtained from the exact bubble to that found 
from various commonly used approximations, such as the Landau-Ginzburg approximation, 
amongst others. We next calculate the activation rate in the presence of finite fermion 
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number. Here we find that beyond a critical baryon number, the activation barriers disappear 
completely. Finally, section 5 contains our conclusions. 



2 The N-flavor Gross-Neveu Model 

We consider the Lagrangian density for the N-flavor Gross-Neveu model in (1 + 1) dimen- 
sions: 

£ = ip a (i^ + m )ip a + ^(4> a 4>a) 2 (1) 

with a = 1 • • • N and a summation convention on a is used. Introducing an auxiliary 
(Hubbard-Stratonovich) scalar field A , the Lagrangian density becomes 

—A 2 N 

C = + V>a + m + A )^a (2) 

For mo = the Gross-Neveu model has a discrete chiral symmetry: 

i>a -> il54>a ; ~> ^a«75 5 A -> -A 

It is this symmetry that will be spontaneously broken (only discrete symmetries can be 
spontaneously broken in 1 + 1 dimensions). The explicit mass term mo introduces an ex- 
plicit chiral symmetry breaking term and will be responsible for a non-degenerate vacuum 
structure. 

We should note that at tree level, the field A has no dynamics; it is truly an auxiliary field. 
However, it acquires dynamical content by dint of quantum effects arising from integrating 
out the fermionsHl. 



2.1 The Gap Equation 

Let us first consider homogeneous configurations A(x) = A and the case of zero temperature. 
In order to understand this case clearly, it is convenient to pass to the Hamiltonian 

r ( A 2 N d 1 

h= y dx {^A7 + ^ HcT3 ^ +aiMo] ^i (3) 

with Mq = m + A being the effective fermion mass, and the cr, are the Pauli matrices. The 
ground state is obtained by filling up the negative energy Dirac sea, so that the ground state 
energy density is given by: 

E Aq 1 f A 



NL 2A 7T Jo 



dkJk 2 + M 2 (4) 
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where L is the spatial length of the system and A an upper momentum cutoff. Dropping a 
term proportional to A 2 (normal ordering) the energy density becomes 

E _ Al Ml Ml ^ /2A_\ (5) 



NL 2A 4vr 2tt \\M ti/ 
while the gap equation (8E/dA = 0) becomes 

A n \\M \J 

Clearly both the energy density and the gap equation need non-trivial renormalizations. We 
find that the proper renormalization conditions are 

M = M R = A R + m R (7) 

f 11, /2A\ 

T-~T ln — (8) 



\r A 7T V K J 

M R m M R ( 2A \ 

? = ? (10) 
Ao Ar 

where k is an arbitrary renormalization scale. The renormalized gap equation now becomes: 

Mr TJIr Mr ( K \ 

"Ar" ~ ~hi = \\Mr~\J (U) 

We see that for m R = the gap equation generates dynamically a scale 

\A g \ = Ke- w/XR (12) 

This solution to the gap equation has the lowest energy (see below), thus spontaneously 
breaking the discrete chiral symmetry. It is convenient to absorb this dynamically generated 
scale as the overall scale in the problem and define all dimensionful quantities in terms of 
this scale. We thus introduce the following dimensionless quantities 

Tiir m r m R 

M = T7H ; m = T7T\ ( 13 ) 

\t-±g\ \L±g\ 

Discarding terms independent of M, the dimensionless energy density becomes 

E M - , . , Mm 

— MM 2 ) - 1] - — (14) 



and the gap equation becomes 



i^-S- (15) 
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Clearly, for m = the solution M = ±1 gives the lowest energy, thus the discrete chiral 
symmetry is indeed spontaneously broken. 

The above energy density for a constant Ar configuration is identified as the fully renor- 
malized effective potential for the auxiliary scalar field. 

At finite temperature T we compute the free energy for the constant configuration and 
thus obtain the finite temperature effective potential. The free energy is given by 

F[A] = -f ln(Tr exp{-H[A]/T)) (16) 

where the trace is over the fermionic degrees of freedom. Discarding irrelevant terms, and 
introducing (along with M, m above) the dimensionless temperature T = T/|A S |, we find 
the dimensionless free energy per unit length to be (see Fig. 1) 



roo 

I(y) = / dxln 
Jo 



1 + e 



-^Jx 2 +y 2 



A small-?/ expansion of I(y) is given by| 



I(y) = ^ + j MvM + 7 - 1/2] - + 0{y") (19) 



where 7 ~ 0.5772 is Euler's constant and ((3) ~ 1.202 is the Riemann zeta function. 
The extrema of the free energy are obtained from the gap equation: 

This equation can either have three solutions (corresponding to a local minimum [metastable 
state] at M — M_ < 0, a maximum, and a global minimum [stable state] at M — M + > of 
the free energy), or just one solution (minimum). The left hand side of ( ^0|) that determines 
the local minimum at M_ is plotted in Fig. 2 for various values of T; we see that solutions 
only exist for m/A# < \/{&k) ~ 0.117, and for a given m only over a certain range of T's 
{e.g. for m = 0, T < e 7 /7r ~ 0.567). 

2.2 The "Bubble" Configuration 

Static topological and non-topological semiclassical solutions in the Gross- Neveu model have 
been analyzed and systematically constructed by Dashen, Hasslacher and Neveu (DHN) HlQ[ 
and Campbell and Liao||ll|| using inverse scattering techniques. These solutions were also 



studied within the context of electron-phonon models of conjugate polymers []T2, |T5| , [TJ], IE 
and we refer the reader to these references for a comprehensive analysis of these solutions. 
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For the case m = (explicit chiral symmetry) there are two degenerate solutions to the 
gap equation. In this case there exist semiclassical topological soliton (kink) solutions and 
non-topological bubble (polaron) solutions [|TU|, |Tl"| , |T2"f . However when m 7^ 0, the effective 
potential (or free energy density) has one local and one global minimum, the degeneracy 
is lifted by the explicit chiral symmetry breaking mass term and no soliton solutions are 
available. However, by taking the results from the degenerate case, it is straightforward to 
find that in the m^O case there still exist bubble solutions for which the "gap" parameter 
A (and therefore M) is spatially varying. These solutions (obtained via inverse scattering 



methods [TT| , |12| ) are given by 

M±(x) = M± — k {tanh [k (x — x + y)] — tanh [k (x — x — y)}} (21) 
where M± are the minima solutions to the gap equation The parameters k and y are 



related by the "integrability condition" flnj, |TI 



tanh[2fc y] = (22) 

Note that we may take k > since the sign of k cancels on both sides of the equation. 
Furthermore, notice that the solution requires that \y\ > 1/\2M±\. This is because the bound 
states cannot be localized in distances smaller than the Compton wavelength of the fermion. 
The position of the center of mass of the bubble is determined by xq. This variable reflects 
the underlying translational invariance of the system. Because of this invariance, the energy 
(or free energy) is independent of this variable. For notational simplicity we will take it to 
be zero in the rest of the analysis. 

The fermion spectrum in presence of this non-topological semiclassical configuration is 



also known exactly |12| , pL3|, |14|, |15(| : N bound states at energies 



±w 0± = ±^M^-k 2 , (23) 



and positive and negative energy continuum states with energies ±Wk± = ±yM| + k 2 and 
phase-shifts 

6(k) = 2 tan" 1 f^j . (24) 



with ko the solution of the integrability equation ( 22f) . The exact fermionic wave-functions 



may be found in reference ||12||. A noteworthy feature of these wave-functions is that the 
bound-state fermionic wave-functions are localized at x = ±y. It is known [|T6|, |i~Tf that in the 



presence of a topological soliton (kink) there are fermionic zero modes whose wave-functions 
are localized at the position of the kink. Since the bubble configuration is reminiscent of 
a kink-antikink pair separated by a distance ~ 2|y|, the localized fermion bound-states are 
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the symmetric and antisymmetric combination of the zero-mode wave-functions, split off in 
energy because of the non-zero overlap of these wave-functions. 

The above configuration, being similar to a kink-antikink pair separated by a distance 
2\y\, can clearly be identified with a "bubble". The value of y (the remaining unspecified 
parameter) may be obtained by extremizing the energy of this configuration as a function of 
this parameter. For the moment we will leave this parameter free and treat it as a "collective 
coordinate" . 

Before proceeding to the numerical analysis of the bubble energy, it proves illuminating 
to understand some features of this configuration. 

Asymptotically this solution reaches M±. For the global minimum M + > and the 
integrability condition (^) implies that y > 0. In this case, in the region — y < x < y, the 
bubble samples a region where the free energy (or energy) density is higher than that of the 
global minimum (for large y, probing the metastable state). Thus, for large y, the energy as 
a function of y will grow linearly (because the difference of the tanh's is roughly constant in 
this region). On the other hand, if the bubble solution has its asymptotics in the metastable 
minimum for which M_ < 0, then y < and the bubble configuration is probing a region 
(y < x < —y) in which the free energy (or energy) density is lower, thus gaining volume 
energy. In this case for large y the energy of the bubble configuration will diminish (linearly 
for large \y\). Certainly for small y the spatial variations of the bubble configuration will 
raise the energy, and since for large y the energy will diminish, there has to be a maximum 
of the energy as a function of y. This maximum, the critical bubble, thus corresponds to 
a "sphaleron" in the sense of Manton and Samols|18|. The maximum of the energy as a 
function of the collective coordinate y corresponds to a saddle point in functional space. 
The unstable coordinate corresponds to Sy = y — y* where y* corresponds to the maximum 
of the energy. The collective coordinate xq corresponding to the center of mass of the bubble 
is a flat direction in functional space, along which the energy is constant as a consequence 
of translational invariance. This configuration is an exact saddle point of the full effective 
action (functional of A (a;)) obtained after integrating out the fermions 

A 2 

J3 

V 

Finite temperature enters in the trace via the boundary conditions on the fermionic fields. In 



S eff = N J d 2 x<-£p> -iNTr\n[i$ + mo + A ]. (25) 



references ([10], [11], [12|]) the effective action is written in terms of the scattering data (bound 
state energies and normalization and phase shifts) and the bag configuration is found by 
extremizing the full effective action with respect to these parameters. (This is the essence 
of the inverse scattering method). In the large N limit the functional integral is determined 
completely by the saddle points. However, even for finite N, the critical droplet configuration 
is an exact saddle point of the effective action. 
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Hereafter we will choose (making use of dimensional transmutation) the renormalization 
scale k to be the value (analogous to A QCD ) for which X R = 1. In units of the dynami- 
cally generated scale, the renormalized free energy difference between the metastable bubble 
configuration and that of the constant "gap" configuration in the metastable state is 
F 



N 



^Jdx (A 2 (x) - A 
1 + e~ W0 -/ T 



- w - 



«/ 2 dA; dS 2k 

——w k _ + 

7i ak 7r 



- 2Tln 



- 2T 



00 dk d5 



In 



1 



-Wk-/T 



7T dA; 

The phase-shift S(k) is given by (0), the bound state energy w _ is given by 
are the continuum energies. F/N then simplifies to 



(26) 
and Wk- 



F 

N 



7i 



ATimy + 2 - In M 



2Tln (l + 



-WQ-/T 




-w k -/T 



k 2 + k 2 



(27) 



The first term inside the square brackets shows that the energy diminishes for large \y\ 
(recall y < 0), consistent with the argument presented above. This is the classical "volume" 
energy contribution and arises primarily from the first (integral) term in (p6|) . The rest 
of the bracketed expression is the T = quantum contribution from the fermionic degrees 
of freedom, and saturates at a constant value for large \y\. The remaining terms give the 
thermal contribution. 

We want to use the above results to calculate the activation rate between the metastable 
and the true vacua. Our procedure will be as follows. We have an exact solution to the full, 
quantum equations of motion. To calculate the activation rate we can use the path integral 



formalism [[0J; what we do is to saturate the path integral with the bubble configuration 
found above and evaluate the rate semiclassically. This is just the standard procedure that 
is commonly used. Thus, when we say we are calculating the "exact" activation rate, we 
mean that we are using the exact solution to the full equations of motion to perform the 
semiclassical calculation of the rate. In the large N approximation, our calculation would 
become exact in all senses of the word. 

In the large N limit the gaussian fluctuations around the saddle point configuration yield 
a contribution of order (1/iV) and to leading order the prefactor is 1. However the free 
energy F = N(F/N), with F/N given by equation (p7|), and the rate formally vanishes in 
this limit. 

For finite N we should take into account the corrections from the gaussian fluctuations 
around the saddle point configuration, that is the prefactor, this is an extremely difficult 
task in this case. 

Therefore, we only concentrate on the exponential since it gives the leading behavior 
and in this case it contains the quantum corrections associated with the fermionic loops, 
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in particular in the case under consideration, the contribution to the exponential is solely 
arising from quantum corrections, since at tree level the auxiliary field has no dynamics. 

Having clarified the nature of the approximations involved, we now set N — 1. 

Given m and T, then, we can plot F as a function of y (Fig. 3), and the maximum 
is F c c . The superscript "C" indicates that this is the "Correct" method of calculating the 
critical bubble free energy. Up to a prefactor, the (semiclassical) activation rate is then 
r (T) = exp(— /T). For m = 0.02, T C (T) is plotted as a solid curve in Figs. 5-8. 



3 Some Approximation Schemes 

Now we turn to the main point of this work, which is to compare the exact critical free 
energy in the GN model to that obtained by means of various approximations. 

There are a variety of ways to approximate the free energy that enters into the activation 
rate. The class of schemes we consider involve approximating the kinetic term in various 
ways (note that the theory written in terms of the auxiliary field A has no canonical kinetic 
term at tree level) as well as approximating the potential term. Some of these schemes will 
have counterparts in the 3 + 1 dimensional calculations of activation rates in scalar fields, 
while others are more specific to the Gross-Neveu model. 



3.1 The Landau- Ginzburg Approximation 

The BCS theory of superconductivity contains a four-fermi interaction which controls the 
formation of Cooper pairs. A complex auxiliary field <fr can be introduced to rewrite this four 
fermi interaction as a Yukawa type coupling. This field then serves as an order parameter for 
the superconducting state. Its dynamics can be determined by integrating out the fermions 
and constructing the effective action for 0. To understand the nature of the phase transition 
in this case, it is useful to perform a derivative expansion of the effective action near the 



critical temperature. This is the Landau- Ginzburg expansion f20[| . We can perform the same 
procedure for the Gross-Neveu model. 

The Landau-Ginzburg effective free energy is obtained as a consistent expansion in small 
order parameter (M) and small gradients, and is valid near the critical temperature and for 
small explicit symmetry breaking fields (small m). In this approximation (method "L") the 
effective free energy density is written up to terms of order (dM(x)/dx) 2 in the derivative 
expansion and up to order M 4 (x) in the non-derivative (potential) terms. The calculation of 
the lowest gradient term is obtained via a Feynman diagram expansion to one loop using the 
imaginary time formulation of finite temperature field theory. It is carried out by writing 
M(x) = M+v(x) with M the homogeneous configuration, and v (x) (the small departure from 



the homogeneous configuration) taken to be a perturbation. Due to the Yukawa coupling of 
M(x) to the fermions, M will be the mass of the fermions in the loop. 
The coefficient of the gradient term is the bracketed expression in: 

p v(p)v(—p) 



2 



MttM 2 \2ttMT \T)\ y 

where I(y) was defined in eq. fll8|). The diagram shows a fermion loop attached (with coupling 
v) to two (truncated) scalar legs carrying spatial momentum p. The field v(x) is taken to be 
independent of the Matsubara frequencies, as well as slowly varying in space, in keeping with 
the philosophy of the derivative expansion. In the T — > limit the bracketed expression is 
1/(24ttM 2 ), while for T > M it is [using the expansion of eq. ([19])] 7C(3)/(32tt 3 T 2 ). We have 
verified this calculation by means of Chan's method of obtaining the derivative expansion [ETJ]. 

Using the quartic approximation (for T ^> M) to the finite T effective potential of eq. ( |17D 
(see Fig. 4) yields the Landau- Ginzburg free energy density (per flavor): 

This free energy density can be extremized with respect to M(x), and it is found that the 
extremum inhomogeneous configuration is a bubble solution p2[ which is of the same form as 



the true extremal bubble, eq. (|2T|), with k fulfilling exactly the same integrability condition 
as in eq. (|22f) . However, in this case, M_ is the local minimum of the Landau- Ginzburg 
potential, and for this extremal bubble, 

•!4^I h M (30) 

With t = k /\M_\, the total free energy of this inhomogeneous configuration in the Landau- 
Ginzburg approximation is 



f l _ 7C(3) 3 
c ~ 4vr3T 2 



-4 2 1-t 2 n-t 

+ ttt + — ^— In 



3 t 2 t 3 Vl + i 



(31) 



The exponential part of the activation rate, T L (T) = exp(— F^/T), is plotted as a dashed 
curve in Fig. 5, along with the correct result T C {T); from this graph we see that method 
"L" overestimates the rate up until the critical temperature where they both become equal 
to unity (the critical temperatures are slightly different). This feature will recur in all the 
approximation schemes we use here, except for the last (zero T bubble). 

3.2 Landau- Ginzburg With T = Gradients 

In method "G", the same Landau- Ginzburg potential is used, but only the T = gradient 
term is taken, i.e. the coefficient in front of the derivative term is l/(247rM 2 ). The critical 
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bubble free energy is then given by 



if = — , R = J24xM1\ (32) 
R' V V 32tt 3 T 2 v ; 

where Fz 1 was given in eq. fH]). T G (T) = exp(-F c G /T) is plotted dashed curve in Fig. 6. 



This method is closer in spirit to the one used in the standard calculations |19| of finite 
temperature activation rates for scalar field theories with canonical kinetic terms. In these 
calculations, the canonical kinetic term is combined with the finite-T effective potential to 
give the free energy relevant to the activation rate. 

3.3 The Effective Potential Approximation 

In the effective potential approxmation (method "P"), we define 

F pot = Jdx [V cS (M(x)) - Kff(M_)] (33) 

where V c s was defined in eq. ([17]). Among configurations of the form eq. (|21f) [and satisfying 
the integrability condition eq. (^2j)], parameterized by the half-width y, we choose y to 
extremize F pot , and the resulting F pot is our approximation to the critical free energy . 
r p (T) = exp(— F c p /T) is plotted as a dashed curve in Fig. 7. 

Since classically there is no kinetic term in our Lagrangian, eq. (0), this method is 
tantamount to ignoring derivative corrections (both quantum and thermal). In real- world 
calculations, such as that of the free energy of a critical bubble or the mass of a sphaleron, 



derivative corrections are often ignored [E3|]; the finite-T effective potential is combined with 



the classical (canonical) kinetic term, and the resulting approximate action is extremized. 



For example, the sphaleron mass at T = is known to be p4 



AttvB 

M SP = , (5 = 1.5-2.7) (34) 

9 

where v is the minimum of the Higgs potential. The sphaleron mass at finite T is then 
approximated by the same formula but with v — > v(T), the minimum of the Higgs effective 
potential [and also g — > g{T)]. Derivative corrections to the action have been ignored. 

3.4 The Zero-T Bubble 

In the zero-T approximation (method "Z"), the approximation scheme is as follows. We start 
by finding the T = critical bubble, a configuration of the form eq. ([H|) and eq. (p2[), with 



M_ the minimum of zero temperature effective potential, V^t(T = 0), choosing y to extremize 
the T = free energy. Then we plug this value of y into the finite-T formula for the free 
energy, eq. (p7|). T Z (T) = exp(— F c z /T) is plotted as a dashed curve in Fig. 8. 
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This is analogous to the method used for critical bubbles in f25 |. There the bubble 
configuration used was the extremum of the classical, T = action. The (1-loop) finite-T 
free energy of this bubble was then calculated (exactly, including "derivative corrections") 
by numerically evaluating the determinant factor. 



The purpose of ref. |23[ was to measure "derivative corrections" by comparing the Zero- 
T calculation (which includes derivative corrections) to the effective potential calculation 
(which does not). This effort was hampered by the fact that the Zero-T calculation is itself 
an approximation (only good at lower temperatures). In the Gross- Neveu model of this 
paper, we can compare both methods to the true result. 



4 Finite Fermion Number 

If sphaleron configurations are relevant in baryon number violating processes, an important 
question to pose is how the presence of a finite baryon number modifies the field profile and 
the free energy barrier of the sphaleron configuration. Within the model under consideration 
we can provide some partial answers to these questions. In this model we identify baryon 
number with fermion number. In order to consider a finite baryon number density at zero 
temperature we would have to allow for the case in which the N positive energy fermion bound 
states are occupied, and the rest of the baryons to be in the positive energy continuum up 
to a Fermi energy (or chemical potential). Because the phase shifts of the continuum states 
depend on the bubble "collective coordinate" y, to obtain the bubble profile for arbitrary y, 
the Fermi energy will have to be adjusted such as to provide the constant baryon density (a 
constraint on the system). 

At finite temperature we would have to introduce a chemical potential and work in the 
grand canonical ensemble, finally fixing the chemical potential (as a function of temperature 
and bubble size) to give the fixed baryon number density (on average). Both situations are 
extremely hard to implement and at the present time lie beyond our capabilities. 

We will content ourselves with considering the somewhat more restricted scenario in 
which we have a fixed baryon number (B), < B < N and the case of zero temperature. The 
reason that this case is somewhat simpler is that we can accommodate the B baryons in the 
available N positive energy bound states. This corresponds to the ground state configuration 
with a finite (positive) baryon number B < N (for a negative baryon number we would have 
to deplete B negative energy bound states). 

Following the steps leading to ( f27|) at T = we find that the free energy becomes 

F B F° B 

AT = AT + iv"° (35) 
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with F° given by equation (B7I) with T = 0. It is convenient to parametrize 

k = |M_|sin(0) 

uj = |M_|cos(0) (36) 
1 + sin(0)" 



4 sm{9) 
< 9<n/2 



sin(0) 



Since y{6) is a monotonically increasing function of 9 it proves more convenient to locate the 
size of the critical bubble by extremizing the free energy with respect to 9. The extremum 
condition leads to 

Tn . . 9 B . 
tan(0) + - = — (37) 



|M_| w 7T 2N 

As we have analyzed in section two, the metastable minimum is available only for 7rra/|M_ | < 
1. Thus we find that for B ^ there are two non-trivial extrema of the energy functional. 
Clearly one corresponds to a minimum and the other to the maximum, i.e. the critical 
bubble. The free energy barrier between the minimum and the maximum becomes smaller 
as B is increased, finally disappearing at a maximum baryon number given by B max where 
B max is found using 9 max = sin -1 ^1 - ^j^j in eq. (0). 

For B > B max the barrier disappears altogether and there is no longer a bubble solution. 
This is indeed a remarkable result if it persists in three dimensions. The physics of this 
phenomenon is clear. By filling up the allowed positive energy bound states the Pauli pressure 
(a consequence of the Pauli exclusion principle) increases. For the minimum which appears 
at a smaller bubble radius this pressure is larger than that for the maximum, which appears 
for a larger bubble radius. Although extrapolation to three dimensions is not warranted, 
we would expect fermionic bound states in the lowest (J = 1/2) partial wave, because 
in this partial wave the problem is essentially one-dimensional (save for the fact that the 
wave function must satisfy proper boundary conditions at the origin). If the effect of these 
fermionic bound states localized at the wall of the bubble is consistent with the behavior just 
found in this 1 + 1 dimensional model, it would imply that the current bounds on sphaleron 
transitions may be a gross underestimate for the rate. Clearly we cannot conclude that this 
happens in the more realistic three dimensional scenario, but we believe that this observation 
may be worthy of consideration there as previous estimates did not take into account the 
fermionic back-reaction onto the sphaleron configuration and their effect on the activation 
barriers. 
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5 Results and Conclusions 



The amazing thing about the GN model is that one can calculate the exact critical bubble 
configuration of the theory and hence we can deduce the true (modulo our statements about 
the semiclassical evaluation of the relevant path integral) value of the activation rate (or at 
least its exponential part, which typically is the most important part of the rate). 

T(T) = e~ Fc/T calculated by our several methods (one exact and four approximations) 
are plotted in Figs. 5-8 (for m = 0.02 and N = 1). The bubbles M(x) (for m = 0.02, N — 1, 
and T = 0.35) used in these methods are plotted in Fig. 9. 

We see that the Landau- Ginzburg method ( "L" ) and the potential method ( "P" ) converge 
nicely to the correct result ("C") at high temperatures. They both correctly predict T(T C ) = 
1, though their respective critical temperatures vary slightly from the correct T c ; we would 
only expect exact agreement if the phase transition were second order (the limit m — > 0). 
The zero-T method ("Z") is a poor approximation at higher temperatures, but is better at 
lower T's. Curiously, all the approximation schemes tend to overestimate the rate. 

Finally, we have to address the question of what this means for realistic field theories in 
3 + 1 dimensions. Since the GN model was so special in that we could find the exact bubbles, 
etc., one might be tempted to think that our results should only hold for this model. While 
we have no evidence that this is not the case, it should be recalled that the GN model serves 
as a very good testing ground for issues such as asymptotic freedom, which does in fact occur 
in 3 + 1 field theories. It is also the inspiration for Nambu-Jona-Lasinio low energy effective 
models for strong interactions. Thus, our results may have a wider range of applicability to 
theories in higher dimensions. In particular, if the various sphaleron configurations in higher 
dimensional theories can sustain fermionic bound states, then we expect that our discussion 
of the activation rate in the presence of finite fermion number should apply. This would be 
of great importance since the reduction of the barrier would have significant consequences 
for the activation rate. 

At any rate, what our results do show is that the standard methods of calculating the 
activation rate need not be all that close to the exact value; if nothing else, our calculation 
can serve as a warning sign to those who might believe their approximations overly much. 
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V eff 




Fig. 1: V eS (M) for m = .02 and T = {0, .25, .35, .447}. 



LHS 

0.12 r 




Fig. 2: Left-hand side (LHS) of the gap equation at T = {0, .25, .35, .45, .55}. 
Dashed line is m/\n = .02, for which the critical temperature is T c = .447. 
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Fig. 3: Bubble free energy F(y) as a function of the half-width y, for m = .02 
and T = {0, .25, .35}. The maximum (for a given T) is the critical free energy F c . 



V ef f 




Fig. 4: True V e g(M) (solid) and the Landau- Ginzburg 
approximation (dashed), for m = .02 and T = .35. 
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Fig. 5: T(T) = e~ F ^/ T for m = .02: solid = C (Correct), 
dashed = L (Landau-Ginzhurg). 




Fig. 7: T(T) = e~ F ^l T for m = .02: solid = C (Correct), dashed = P (Potential). 
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Fig. 9: Critical bubbles M(x) for m = .02 and T = .35, five methods: 
solid = C (Correct), long-dash = L (Landau- Ginzburg), 

dots = G (Landau-Ginzburg With T = Gradients), 
dot-dash = P (Potential), short-dash = Z (T = bubble). 
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